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ABSTRACT 
Recently, exponential models of refractivity have been 
devised by several investigators to replace the linear, "standard" 
refractive atmosphere. In this study, average values of N are 
evaluated from the data at the mandatory radiosonde levels up 
to 25 mb,and mean coefficients c; satisfying N,,,=\, 2° iii) 
have been evaluated for a group of eight mid-latitude stations, 
and for the entire winter season of 1958-59. The coefficients 


c., plotted as functions of the height of the i-th layer- 


1? 
center, possess certain characteristic features in common at all 
stations. 

The writer wishes to express hisS appreciation for the 
assistance and enoouragement given him by Professor Frank L. 


Martin of the U. S. Naval Postgraduate School in this investiga- 


tion. 
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1. Three month mean-values of Nj at indicated pressure 
levels, and stations employed in this study 
(in N - units) 





RISE OF SYMBOLS USED 
Index of refraction 
Refractivity 
Temperature in degrees Kelvin 
Pressure in millibars 
vapor pressure in millibars 
Height in meters 
Autoconvective lapse rate 
Atmospheric lapse rate 
Curvature of a ray 


Ray-direction relative to the spherically horizontal 
level 


Radius of the earth 


Empirical constant, affording best fit of an 
exponential formula 


Mean value of c satisfying /V.,, = WN, a 7) 
where Z.,,- 2; is the thickness of the i-th standard 
layer 





Ll. Introduction 

The refractive index n is defined as the ratio of the 
speed of an electromagnetic wave in a vacuum to the wave speed 
in the medium concerned. However, this particular form of the 
refractive index is not Vex oomPpericaut in Cheecereebat ion of 
ray-paths; hence the modified refractive index, N is normally 
used. 

The refractivity N, which is equal to (n-l) x 10°, may 
be obtained from radiosonde observations using the equation 

fife A.(p 182 ) me pepwletekO © 
i T (1) 

where T is the temperature in degrees Kelvin, p is the 
pressure in millibars, e is the vapor pressure in mbs, and 
A and 8 are constants. 

Using the Smith & Weintraub [4] values of these 


constants, Equation (1) may be written: 


P m2 5G 
N+ 26h , 3232862 @) 


A "standard" distribution of N in the vertical may be 
obtained by averaging over a sufficient number of individual 
soundings. Bean and Thayer [1] have pieecsted that its” 
curve can be fitted by an exponential formula. The model 
which has recently been adopted by the National Bureau of 


Standards is of form: 


Nf > [we ena (3) 


where N is the refractive index at any altitude, 





N, is the mean surface value of the refracti: | Pie Ae’, 


zis the height above the surface, 


c is an empirical constant, affording "best fit" in 


Equation (3). 

Equation (3) is used only with the mean values of N 
averaged over a long enough period to eliminate all transitory 
trapping layers. The study by Bean and Thayer iy also suggested 
further study by obtaining individual values of Ng and ¢ for 
varying climatic regions. This was the objective of the 
first approach of this author to the problem. Later, however, 
it was decided to make a finer-scale determination of the 
parameter c by determining a mean-value c; for each layer 
between mandatory levels. Thus a working equation of form 


eer ergo re 


// (eee TAS ay Ail ae AP (4) 


: = 
was employed. Here cj is temporarily assumed to be constant 
Ee tne i-th layer, which extends from z; to zj4] and Nj is 
the value of N at z = zj Note finally that cj may be 


obtained by differentiation of (4) as 


- . 1 BN 
a Tae 5) , Gy 





2. Data Processing 

The stations utilized in this study were selected to 
represent a wide variety of climatic regions. The cj~values 
were determined from the three-month means of pressure, tempera- 
ture, and moisture content, for all the mandatory levels at 
these stations during the winter of 1958-59. These data were 
obtained from the radiosonde checked data found in the 
Northern Hemisphere Data Tabulations Daily Bulletins [2] : 

The use of data obtained from radiosondes introduced some 
uncertainty of the relative humidity when 'motorboating" was 
reported. However, motorboating usually occurred near 400 mb, 
and there the error introduced by omitting "motorboating" 
reports from the moisture-.veraging process was found to be 
<€ 0.1 mb. The twice~daily observations taken at 0000 and 
1200 Greenwich times were averaged. This served to filter 
out the diurnal variation, and to make the stations, some of 
which were separated by as much as six hours time-difference, 
more compatible. Two methods were used in obtaining the 
values utilized for the determination of the surface refractive 
index data. When the surface report was part of the published 
radiosonde, this was the preferred data-source. When this 
radiosonde level was not available, it was neeeseaee to employ 
the surface observation for the station. However the surface 
observation was recorded only once each day, at 1200Z, and 


this would result in a comparison of a onceedaily observation 





egainst a twice daily one when going from the surface to 850 
millibars. As a test, the surface three-month mean N-value for 
Washington was computed using both methods and a difference of 
only one N-unit in 300 was noted. Since Washington had one of 
the largest surface diurnal variations of the stations studicd, 
this indicated that the once-a-day surface observation reports 
could be substituted for the two-a-day radiosonde surface 
reports without serious error. 

For the actual determination of mean values, the standard 
levels from the surface up to 25 mbs were used. The twice-daily 
observations at these levels were totaled and then divided by 
the total number of entries. All the calculations were done on 
either tape adding machines or hand calculators, A maximum of 
180 observations were used for each levé1, provided ‘that-un dvta 
were niseisg. The mean vapor pressure was obtained by determining 
the vapor pressure from each observation using che temperature 
and relative humidity data, and finding the mean of these. For 
the majority of the stations, moisture was recorded only up to 
the 400-mb level, at which level the vapor pressure was found to 
be less than 1 mb. The only station that had moisture data for 
the 300-mb level showed the mean vapor pressure to be less than 
one-tenth of a mb. This showed that very little accuracy would 
be lost by not including moisture considerations above 400 mbs. 

After the three-month mean values of pressure, temperature, 
and vapor pressure were obtained for the mandatory levels of 


each station, the three-month mean N; values were determined 





for these levels using Equation (2). These values and the 
appropriate pressure levels are given in Table l. 
The cj -values for each layer were obtained using the 


refractivity equation in the form: 


ie eee 
Pla es ZAG 


where Z;,} "Zi 18 the average-value of the thickness in the 
standard layers employed in this study. 

The values of cy for each station were then plotted on two 
cycle semi-logarithmic paper versus the midpoint-~height of each 
layer as ordinate. These midpoints were taken from the U. S. 
Standard Atmosphere elevations for the stated pressure. A 
discussion of these results and of some applications are given 


ineethe toOllowang sections. 






a & = =_ —/_ 
fa —_———ae 

= 

——_ —_! — 


> —-_ = 


= a 











= SS 


PX , Adak Denver Ship P Tatoosh Kodiak Washington Ship N Lake Charles 


Pressure 3 Meters* | 1611 Meters] 6 Meters 31.Meters 8 Meters 88 Meters 6 Meters 5 Meters 
Surface 313 254 314 323 Ss aio 328 345 
(1003.1mb) |(831.7mb) |(992.6mb) |(1014.8mb) | (1002.8mb) |(1009.7mb)} (1021.mb) | (1020.0mb) 
850 262 ee 261 266 262 261 268 272 
700 216 212 214 217 — 7 217 oie 216 
500 161 *) Bis7®. 158 * 159 3 i161 158 156 154 
400 ~ (3G 13) 131 138 | wee 130 128 i 128 
300 | 105 103 103 104 = 105 ja 102. 100 101 
250 .- 100 88 87 89 _ 88 88 87. 8H 
200 71 72 70 > 2 > 70 om 72 Pls 
150 525 54 |g (52 53 52 54 ‘Pe Ws: 
100 35 36 36 er a | 36 w |: ~ 38 
50 ° > ay 18 L7 os 2) a 7 ee cere. ie. | 
"25 9 9 - 9 9 9 Fe: Om 


~ s - 


* Surface height is listed under each station; station pressure is enclosed in parentheses. 
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Table 1. Three month mean-values of N; at indicated pressure -levels, and stations employed in this 
study(in N-units): | : | 


> 





3. Discussion ard Results 
(a) Typical Model 

By Equation (5), cj is a good indication of the vertical 
distribution of N. An increasing value of cj indicates a 
larger rate of decrease of N with height. Most stations 


studied showed two main values of c- in the vertical. A 


i 
nearly constant value was indicated for the troposphere, and a 
somewhat larger value shown for the stratosphere. These 
values were separated at most stations by a characteristic 
feature, wrich will be referred to as a dual discontinuity. 
The dual discontinuity is evidently associated with the tropo- 
pause; it can be described as having a minimum value of cj 

and immediately above it a pronounced increase to a maximum 
value. A model of cj vs height which depicts schematically 
the association between the dual discontinuity and the 

related temperature sounding is shown in Figure 1, which may 
be regarded as typical of most mid-latitude stations. 


In order to understand the association depicted in 


Figure 1, consider N as given by Equation (2) with e=0. Then 


S77 Glee 
N = 776(2) ” 
Partial differentiation of N in (7) with respect to z, 


and use of the hydrostatic equation, leads to 


ON 2. 2 77 
oo ee ae) 
where &%, is the autoconvective lapse rate, and Y is the 
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atmospheric lapse rate. By use of Equation (5), it follows 
that 


c,e x (4-5) (8) 


Hence c; is relatively large in a thermally stable layer 


1 
and relatively small in a thermally unstable layer. Thus if 
we disregard moisture, any minimum value of cj is associated 
with a relatively steep lapse in the layer, and a maximum 
value with a relatively stable lapse. The discussion of the 
individual stations, which follows, will point out some 
deviations from the typical model, Figure l. 

One can also qualitatively estimate what the effect of 
moisture lapse would be on cj. Empirically one knows that a 
large negative value of AN /dzZ may be associated with a 
strong hydrolapse and this would cause c; to be large 
{ Equation (5)] - On the other hand, if the temperature 
lapse is not abnormal, cj may be a minimum as a result of a 
small lapse (or possibly an increase with height) of water 
vapor. Moisture effects may be associated with some of the 
deviations of the individual c; graphs from the ideal one 
of Figure l. 

(b) Adak 

The tropospheric distribution of c; vs height, 

Figure 2, showed a gradual decrease of cj with increasing 


height up to 9 km. Above this level, a pronounced dual dis- 


continuity was found, indicating the presence of a tropopause 
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layer. A second dual discontinuity of lesser magnitude was 
noted above the first. This was taken as an indication of a 
multiple tropopause, which is in keeping with the model of 
Figure l. 

(c) Denver 

The graph for Denver, Figure 3, shows a smaller value of 
c; in the troposphere than in the stratosphere. There is a 
sharp increase in cj to a maximum value near 13 km, with no 
definite dual-discontinuity minimum below. 

(d) Ship "P'', Tatoosh Island, Kodiak, Washipecoue Dac. 

The graphical results for these stations are shown in 
Figures 4, 5, 6, 7, respectively. These graphs have been 
grouped together because they all have features closely re- 
sembling the idealized model of Figure 1. The tropospheric 
cj; decreases slightly with elevation, and is followed by a 
faster rate of decrease ending in the region of 8-10 km. 
The stratospheric cj; then conforms rather closely to the 
schematic model. Kodiak shows a lesser tendency to have the 
pronounced minimum cj as part of its dual discontinuity. 
This might well be a characteristic of its higher latitude, 
together with a greater variability of tropopause height. 
fe) = ship 'N" 

The graphical results for this station, Figure 8, shows 
a greater decrease of tropospheric cj with elevation than is 
shown in the idealized model. This is due primarily to the 


greater moisture content, and consequently greater N-value, in 
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low levels at this subtropical maritime station. Thus the low- 
level values of cj are larger than those for most other stations, 


although the minimum just below the tropopause has nearly the 


same value (about 11.5 x 10% per meter). 
({) Lake Charles 


The graph for Taj: Charies, Fagurer9, shows Chestxepo- 


f 
Ni 


spheric cj’ -valuc decreasing to a minimum which extends from 


+7 


6 km to 10 km. Above this layer, cj increases -0 {gay 
6 


constant maximum value in the stratosphere. The <3 mire tant 


Ky 


c; in the range 6-10 km is a characteristic of variable tropo- 


o/ 
rae 


pas . neight, during the observation period (c,f. Kodiak, 


Beure 6)! 


(g) Composite of all stations 


sd 
’ 
4 


The composite graph of all stations, Figure 10, conforms 
closely to the idealized model. The cj -value decreases 

slowly with height in the troposphere, and is separated from 
the larger stratospheric czy -values by the dual discontinuity. 
A secondary tropopause is also indicated by the less pronounced 


upper dual discontinuity in the neighborhood of 15 kn, 


within the stratosphere. 


20 





- 
7 
a 
; 
7 


— = 











a~ 

- —_ allie i 
L : aa « fe | 
on —- « : ome — — 


i. | 








—_— =) oe bal = - 
—e— —— + es - 

— te Fn — sh 

a — — ee << — 

a? © =e ee 9 ata” — = 

-— ee lee — 








[oe ..49 | <-m» «€ 
ie et 6 ee ee 
———e a «| =_ee eo we =e 

Ol eee Pee thee! Oke © 
Cert) Fe ee ee ——ee ee 
—> weet «O18. Cc. 
A Ore eee ee oe ee lee 
~=—43 eS «& 





4. Other A plications and Conclusions 

By means of Equation (6) a method of solving for Cj 
between mandatory radiosonde levels was established, ‘ie 
resulting value of cj; was there the mean value for the layer 
of thickness AZ;,*2;,,- Z, +» However c; turns out to be 
somewhat variable with height. Thus in Figures 2 through 
9, all graphs of c; vs height have been drawn as continuous 
functions. This means that Equation (4) may be reconsidered 
from the viewpoint that cj be variable with height in any 

y 

layer. In this case, however, the mean value Ci between 
standard levels shall be the same as that computed by 
Equation (6). 


In the field of radar-propagation, one frequently must 


deal with the ray-curvature K=2%%s5 given by 


K = 10°° aoa +018) cm 2 (9) 
where / is the ray-direction relative to the spherically 
horizontal level iar Up to several years ago, it was 
customary to use as a standard ON/az,4 value characteristic 
of the low troposphere, namely an fyzF - si Cx where Q 
is the radius of earth. We know now from the work of Bean 
and Thayer and others (including this work) that (AN Zo )avg 
is negative but decreases in magnitude. Thus, by Equation 
(5) with c; and N both decreasing with height in the tropo- 
sphere, Jama] would be relatively small just below the 


tropopause, 


Consider still another refinement in Equation (4), 
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name ly, Chat c; is empirically iven es a function oP iherent 
of the type shown in Figures 2 through 9, inclusive, for the 


localities considered. Then logarithmic differentiation of Equation 





(4) gives | gy ee 
7 
22 Nos - Cerne Tots _ 6 Ce 
all 74 BS | ie A Ley io (2-2) (10) 


This expression for @n/qgz has a correction tac.us, as 
compared with Equation (5), given by the second term on the 
right side of (10). Actually this term can be readily com- 
puted in finite differences from graphs of Cc; 5S HOG eee, ) 
(for example, see Figures 2 through 9). The correction 
factor may take on some importance in radar-propagation in 
the tropopause layer, in view of the large increase Gf cy 
wieheewevation cecurring there, 

Finally it should be noted that if an equation of the 
Bean and Thayer type, Equation (3), is desired. It may be 


shown that c is a weighted-mean of the individual C;, Or 


PA ie = CC, AZ 9 ue Orne (11) 
Equation (11) makes it apparent that the use of a mean 
c for the troposphere is reasonable. However, it suggests 
that one should use a different weighted-mean for the lower 
stratosphere than for the troposphere. In Equation (11), 
nets one Aerie above the surface, so that the averaging 
must extend down to the earth's surface. Equation (11) may 


be readily solved fo. c using any graph:efsthe type shown in 


Ze 










Oo en — ee 
_ 







= < ae — a, 
= 2.2 2s eepP 
eh eee ee 








Figures 2 through 9, by selecting the c;-values and weighting 
them by the indicated thicknesses AZ;. However, it should be 
remembered that in most cases only half of the lowest layer, 
OZ, , is shown in these figures, For the stations studied 
in this paper, the Ng -values appropriate to the Bean and 
Thayer exponential form, Equation (3), are listed in the first 


row of Table l. 
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